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SUMMARY 


Results are presented of a flight investigation conducted to survey 
the flow field generated by airplanes flying at supersonic speeds. The 
pressure signatures of an F-100, an F-lOl, and a B -58 airplane, repre- 
senting widely varying configurations, at distances from 120 to L 25 feet 
from the generating aircraft and at Mach numbers from 1.2 to 1.8 are 
shown . 

Calculations were made by using Whitham's method and were compared 
with the experimental results . The procedure used for calculating 
Whitham's F(y) function is given in the appendix. Good agreement was 
obtained for the bow shock-wave strength, and overall characteristics of 
the flow field were predicted; however, all details of the flow field 
were not accurately predicted. This was especially true at the higher 
Mach numbers. Since the area distribution used in the calculations was 
for a Mach 1 equivalent body of revolution, experiments were made to 
obtain the Mach line area distribution for one of the configurations. 
When this area distribution was used in the calculations, the results 
were improved. 


INTRODUCTION 


The flow field generated by an airplane flying at supersonic speeds 
can be separated into three regions: local, far, and near fields. The 

local flow field extends from the body surface to approximately the wing 
tips. The flow in this region is important in determining the inter- 
ference effects of one airplane component on another component. The far 
field is that region where all the intermediate shocks have joined and 
only the bow shock and rear shock remain. This flow field is especially 
important in estimating the noise on the ground produced by supersonic 
airplanes flying at high altitudes. The near field is the region beyond 
the wing tips of the airplane and including the intermediate shocks pro- 
duced by the canopy and wings. Determination of the flow pattern in this 
area is necessary in assessing the danger involved in close-formation 
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flying, aerial-refueling operations, and other situations involving 
close approaches or passes. In reference 1 it was shown that, at 
supersonic speeds, appreciable motions can be induced on an airplane 
flying in proximity to another airplane. In addition, when an airplane 
is flown at supersonic speeds at extremely low altitudes, it is neces- 
sary to take into account the entire shock pattern to predict ground 
pressures and damage to ground establishments. 

Much theoretical work has been done to predict these flow fields. 
The method of characteristics (ref. 2) has been used with considerable 
success in predicting local-flow characteristics. The method developed 
by Whitham (ref. 3) is applicable to both the near field and the far 
field. This method has been used previously (refs. 4 to 6 ) to predict 
the far field and has been shown to give excellent agreement with 
experimental data. 

This paper is concerned with the flow in the near field and pre- 
sents both experimental and theoretical data. Three airplanes (fig. 1), 
an F-100, an F-104, and a B- 58 , were used in these tests to obtain the 
experimental flight data. Pressure measurements were made at several 
distances from each airplane to determine the effects of body shape and 
distance away from the airplane on the shape and intensity of the shock 
pattern. In addition, Whitham 1 s method was used to determine if these 
results could be accurately predicted. 


SYMBOLS 


F(y) = — s "(t)dt 

2 « Jo jTTt 

i,m,n integers 
k = 2 _1 / 2 (7 + 

l length of airplane 

length of flow field 

M Mach number 

Pq static pressure, lb/sq ft 

Ap incremental pressure, lb/sq ft 
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S cross-sectional area 

t variable 

x longitudinal coordinate 

y,r cylindrical coordinates (y, measured along body axis; 

r, radial coordinate perpendicular to y) 

3 = n/m 2 - 1 

7 ratio of specific heats for air (l.4) 

A prime (') is used to denote the first derivative of a function; 
a double prime (") denotes second derivatives. 


DESCRIPTION OF TESTS 


The flight-test procedure used in this investigation to obtain the 
experimental data on supersonic flow fields is illustrated in figure 2. 

The flow field to be measured was generated by the lead airplane, 
which was flown at a preselected Mach number. The passing airplane was 
flown through this flow field at a constant passing rate, and variations 
in static pressure were measured by a sensitive pressure transducer 
mounted on the nose boom of the passing airplane. To determine the 
lateral-separation distance, a motion-picture camera was mounted on the 
passing airplane to photograph the lead airplane. 

Data were obtained on the F-100 and F-104 airplanes at a Mach number 
of 1.2 and on the B -58 airplane at Mach numbers from 1.3 to 1.8. Lateral- 
separation distances were from 120 to 425 feet. 


RESULTS AND DISCUSSION 


A simplified illustration of a supersonic flow field is given in 
figure 3 . The effects of airplane configuration, separation distance, 
and Mach number on this flow field can be seen by examining figures 4 to 
6 . Because of the experimental procedure used, the absolute lengths of 
the flow fields were not determined; therefore, the measured pressures 
shown in these figures are presented in terms of nondimens ional distance 
between the front and rear shocks. 
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The effect of airplane configuration can best be seen from fig- 
ures 4(a) and 5 (a), since data for both figures were obtained at a Mach 
number of 1.2 and a separation distance of l60 feet. The fineness ratios 
of the two airplanes are approximately the same and it can be seen that 
the peak bow-wave pressures are also approximately equal. Configuration 
differences are reflected in the location and strength of the shocks 
following the bow shock. 

The effects of separation distance are apparent in figures 4(a), 

(b) , and (c) and in figures 5(a) and (b) . When the separation distance 
is increased, the shocks decrease in strength and move outward until, 
eventually, all the shocks have joined with the bow and rear shocks to 
form an N- shaped wave. This gradual formation of the N-wave is espe- 
cially evident in figures 4(a), (b), and (c); in figure 4(c), all shocks 
except one intermediate shock have joined with the bow and rear shocks. 

Increasing the Mach number from l.J to 1.8 (figs. 6 (a), (b), and 

(c) ) has little effect on the flow field. Theoretically, the bow shock- 
wave pressure increases slightly with increasing Mach number; however, 
at these Mach numbers and separation distances this is not obvious, 
since the differences in bow-wave pressures shown in these figures are 
approximately those predicted for such differences in lateral separation. 

To calculate the flow-field pressures, the airplane was replaced by 
its Mach 1 equivalent body of revolution using the normal cross-sectional- 
area distribution. The area distributions for the airplanes are shown in 
figures 7(a), (b), and (c). The flow field was then calculated for this 
equivalent body by using Whitham's method as outlined in the appendix. 

This calculated flow field is shown in figures 4 to 6 by the solid lines. 

Examination of figures 4 to 6 shows that the bow-wave pressure is 
predicted fairly well for all three airplanes. However, agreement is 
better for the F-104 and B -58 than for the F-100, possibly because the 
F-100 has an open nose and the effective equivalent body of revolution 
is not estimated as well. The general trends of the flow field are seen 
in the calculations; however, the additional shocks that originate from 
airplane components such as the canopy and wings are not accurately pre- 
dicted either in magnitude or location. In particular, the B -58 flow- 
field calculation (figs. 6 (a) to 6 (c)) shows a large negative peak in 
pressure ratio near the back of the flow field, whereas the experimental 
data show the negative peak nearly in the center of the flow field. 

A possible explanation for this discrepancy in the B -58 flow field 
is the unusually large percentage of cross-sectional area away from the 
airplane centerline, since four nacelles are mounted on the B -58 wing. 

For the F-100 and F-104 airplanes, however, the cross-sectional area is 
more nearly centered on the airplane axis. Also, the Mach numbers at 
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which the B -^8 was flown were further from the Mach 1 assumption made in 
determining the equivalent body of revolution. 

Another calculation of the B -58 flow field was made, therefore, by 
using an equivalent body of revolution based on the area intercepted by 
a series of parallel planes swept at the Mach angle and perpendicular to 
the longitudinal plane, as shown in figure 7 (c) compared with the normal 
distribution. The results of thi^ flow-field calculation gave much better 
agreement (fig. 8 ) with empirical data than was obtained by the previous 
method. This concept was also used in reference 7 for calculating the 
wing-body drag at supersonic speeds. This slant cross-sectional area was 
obtained experimentally by immersing a model in water and measuring the 
volume displaced. These measurements were taken every 0.2 inch, which 
corresponds to a body length of y/l = 0.0114. Special attention was 
given to the problem of eliminating any errors caused by the surface 
tension of the liquid. Previous work had indicated that when a small 
model is used these errors are significant if the incremental volume 
displacement is measured by any method requiring liquid flow. To elimi- 
nate this possible source of error, displacement rods were used which 
were extracted slowly as the model was immersed, thereby keeping the 
liquid level constant. When the results of these measurements were 
plotted, the maximum scatter was 7 percent of the maximum cross-sectional 
area. 


CONCLUSIONS 


The results of flight measurements and calculations of supersonic 
flow fields indicate that: 

1. The strength of the bow shock wave and the overall characteris- 
tics of the flow field can be estimated fairly well by using simple 
theory; however, the location and magnitudes of all the intermediate 
shocks are not accurately predicted. 

2. The normal -area distribution used in the theory may not be an 
adequate representation of an equivalent body of revolution for esti- 
mating the entire flow field, especially at Mach numbers much greater 
than 1. Using an equivalent body of revolution based on the area inter- 
cepted by parallel planes swept at the Mach angle greatly improved the 
results of the calculations . 


Flight Research Center, 

National Aeronautics and Space Administration, 
Edwards, Calif., August 22, i 960 . 
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APPENDIX 

PROCEDURE USED FOR CALCULATING FLOW FIELD 


The theory developed by Whitham is given in. detail in reference 3* 
A numerical solution of the F(y) integral appearing in the theory is 
given in this appendix, along with a brief summary of the method used 
to calculate the flow field. 

From reference 3 


F(y) 

p 0 ^23 r 1 / 2 

(1) 

pr - kF^r 1 / 2 + y 

(2) 


Where 


k = 2 _1 / 2 (r + l)M J+ p -5 / 2 


and F(y) is a function determined by the body shape and for a smooth 
slender body is 


F(y) = P S "^ - dt (3) 

231 Jo ^y - t 

If no analytical expression is known for the area distribution, 

F(y) must be integrated numerically. Equation ( 3 ) is a form of Abel's 
integral equation and, because of the singularity of the integrand at 
the upper limit, the transform 


S'(t) = 2 P gfr.fo 
Jo v/t - y 


00 


is used to solve numerically for F(y) in the following manner. 
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Assume F'(y) is a constant over a small interval Ay. Then F(y) 
can be written 


F(y) = F(y) ± + F'(y) i+1 (y _ y.) (5) 


over the interval 


yt = y = yi+i 


where 


i = 0, 1, 2 , , . . m 


yo = F(y ) 0 = 0 




1=0 Jyi 


yn+i " y 


The integral in equation (6) can easily be integrated to give 

^'(yn+i) = II - 2 ( F i - F, i+iyi)[(y n+ i - yi+i ) l/2 - (y n +i - yi) L / ; 


+ 2F '. +1 - fe 1 — 71+1 - y n+ i(y n+ i - y i+ i ) l/2 


~ n+1 - + y n+ i(yn+i - yi) 1/2 
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The only unknown in equation ( 7 ) is F' n+ ^. When this is found, the 

next value for F(y) is known. By this step process the entire F(y)- 
curve can be constructed. 

Once the F(y) -curve is known, it is then simple to find the shocks. 
Whenever F'(y) >0 a shock will appear. To find the bow shock, an 
integrator can be used and, by trial and error, values of y and F(y) 
can be found for which 


— kr 1 /%^( y) = j F(t)dt (8) 

2 Jo 

These values are then used in equations (l) and (2) to determine the 
position and strength of the how shock. 

To find the other shocks, lines with slopes which are given by 

F(y 2 ) - F( yi ) __ _j_ 
y 2 - n kr 1 / 2 

are drawn through the F(y) -curve such that "lobes" cut off on each side 
of the curve are equal in area (see following sketch). These, then, are 
the values of F(y) and y to be used in equations (l) and (2). 
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(a) M = 1.3, r = 150 feet. 
Figure 6.- B-58 flow field. 




(b) M = 1.6, r = 180 feet. 


Figure 6.- Continued 
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